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Abstract 

Summarized  is  a  method  to  count  the  number  of  invariants  in  a  lattice- 
gas  automata.  A  simple  calculation  proves  that  a  deterministic  FHP 
lattice-gas  with  a  definite  chirality,  possesses  only  three  invariants:  the 
total  mass  and  the  two  components  of  momentum.  This  model  in  some¬ 
what  simpler  than  FHP  and  is  time  reversible. 


1  Introduction 

The  discovery  that  a  very  simple  cellular  automaton  model,  the  Frisch,  Hass- 
lacher,  and  Pomeau  (FHP)  model  [1],  using  only  six  bits  of  memory  at  each 
point  in  space  reproduces  hydrodynamics  in  its  macroscopic  limit  has  stim¬ 
ulated  many  investigations  of  the  new  domain  of  discrete  kinetic  equations. 
This  domain,  broadly  termed  lattice-gas  dynamics,  has  to  date  expanded  its 
scope  beyond  incompressible  hydrodynamics  to  thermodynamics,  multiphase 
phenomena,  magnetohydrodynamics,  reaction-diffusion  systems,  and  seems  to 
be  continuing  its  growth.  Given  the  conceptual  importance  of  a  simple  lattice- 
gas  model  with  its  ability  to  simulate  so  many,  the  question  arises  if  their  exists 
a  lattice-gas  cellular  automaton  model  even  simpler  than  FHP  that  yet  has  no 
more  than  three  invariants:  the  total  mass  and  the  two  components  of  momen¬ 
tum.  A  very  simple  calculation,  similar  to  that  previous  done  [2],  shows  there 
indeed  exists  a  simpler  model. 

It  is  well  known  that  symmetric  3-body  collisions  must  be  included  along 
with  the  FHP  model’s  even  and  odd  chirality  2-body  collisions  to  achieve  the 
correct  macroscopic  limit1.  However,  the  fact  that  both  even  and  odd  2-body 
collision  possibilities  are  present,  dictates  when  a  2-body  collision  occurs,  a  coin 
toss  must  also  occur  to  determine  the  outgoing  state.  This  coin  toss  makes  the 
FHP  model  irreversible.  Clearly  if  the  model  is  constrained  to  have  a  definite 
chirality,  say  for  example  2-body  collisions  generated  by  ^  rotations  are  retained 
while  ^  collisions  are  discarded,  then  the  lattice-gas  would  be  strictly  reversible. 
But  does  this  additional  constraint  of  a  fixed  chirality  engender  any  spurious 
invariants?  The  answer  is  no,  see  the  results  given  in  §3.2  for  the  proof. 

1  Antisymmetric  3-body  collisions,  2-body  collision  with  a  spectator  particle,  could  also  be 
included  but  this  is  not  necessary. 
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2  Some  Preliminaries 

Symbols  used 


Spatial  Unit  : 

l 

Temporal  Unit  : 

T 

Particle  Speed  : 

l 

C  —  — 

T 

Momentum  Directions  : 

B 

Lattice  Vectors  : 

ea 

a  = 

1,2,  ...,B 

Distribution  Function  : 

fa 

Collision  Operator  : 

Oa 

Jacobian  Matrix  : 

dab 

It  should  be  understood  that  whenever  the  single  particle  distribution  function 
is  written,  its  subscripted  index  is  taken  modulo  B 

fa+b  f  mod  b(^+6)-  (1) 

The  lattice  Boltzmann  equation  is 

dtfa  +  ceaidifa  =  tta-  (2) 

A  general  collision  operator  is  constructed  as  follows 

=  5>Qa({CJ),  (3) 

{Cd 


where  {£;}  is  a  set  of  occupied  particle  states  and  a  =  ±1  is  a  scalar  coefficient 
and  where  each  term  in  the  sum  is  written  in  factorized  form  as 


Qa (^1  >  •  *  •  ?  ^k) 


fa+il  fa+ik 

1  —  .fa+ii  1  —  fa+zk 


(! -/»+;)• 
j= 1 


We  expand  the  distribution  function  about  its  equilibrium  value,  fe 

fa  =  r  +  6  fa 

so  that,  to  first  order,  we  have 

Va{n  =  YJffjh- 


(4) 

(5) 

(6) 


The  l.h.s.  of  (6)  must  vanish,  since  the  particle  distribution  is  non-changing 
under  equilibrium  conditions.  The  eigenvalues  of  the  Jacobian  of  the  collision 
operator, 


Jab 


Ottg 

dfb ’ 


(7) 
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can  be  calculated  and  the  number  of  these  that  vanish  must  equal  the  number 
of  invariant  quantities  in  the  lattice-gas  dynamics.  Because  of  the  finite-point 
group  symmetry  of  the  spatial  lattice,  the  Jacobian  matrix  will  be  circulant  and 
therefore  its  elements  can  be  specified  by  the  difference  of  their  indices,  Jab  = 
Ja-b-  This  property  of  the  Jacobian  simplifies  the  solution  of  the  eigenvalue 
equation 

J2ja-btf  =  \k&  (8) 

b 

where  k  =  1 , . . .  ,B.  Let  us  make  the  ansatz  that  the  eigenvectors  have  the 
following  form 

£  =  ^iak/B  (9) 

Then  inserting  (9)  into  (8)  and  taking  m  =  a  —  b,  gives 

\k  =  ^  Jme2nimk/B.  (10) 

m 

3  Triangular  Lattice:  B=6 

3.1  Eigenvectors 

Using  (9),  the  eigenvectors  of  the  Jacobian  matrix  are 


6  = 

(i,i,  i,i,  i,i) 

(11) 

a  = 

(e,e*,-l,e,e*,l) 

(12) 

a  = 

(e*,e,l,e*,e,l) 

(13) 

a  = 

(-1,1, -1,1, -1,1) 

(14) 

a  = 

(e,  e* ,  1,  e,  e*,  1) 

(15) 

a  = 

(e*,c,-l,e*,e,l), 

(16) 

where  e  =  e*  3 . 

3.2  Definite  Chirality  Models 

Even  chirality: 

Ua  =  Qa(l, 4)  -  Qa(0, 3)  +  Qa{  1, 3, 5)  -  Qa{ 0, 2, 4)  (17) 

fio  =  -  (/o  (1  -  /i)  (i  -  h)  h  (i  -  h)  (1  -  h))  + 

(1  -  /o)  h  (1  -  /a)  (1  -  h)  U  (1  -  h)  - 

fo  (1  -  /i)  /a  (1  -  /a)  U  (1  -  h)  + 

(1  -  f0)  fi  (1  -  h)  h  (1  -  U)  h 

j  =  drC[(-1  +  /)2/2,-((-1  +  /)2/2),(-1  +  /)2/3 

(1  -  2  /)  (-1  +  /)2  /2,  (-1  +  /)2  /2  (-1  +  2  /) ,  -  ((-1  +  / )2  /3)  ] 
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A0  =  0 

Ai  =  0 

A2  =  2e(l  +  e)  (-1  +  /)3/2 

A3  =  -6  (-1  +  /)2/3 
A4  =  2e(l  +  e)  (l-/)3/2 

A5  =  0 

There  are  only  three  zero  eigenvalues,  so  this  even  chirality  model  is  sufficient. 
Odd  chirality: 

fia  =  Qa( 2, 5)  -  Qa( 0, 3)  +  Qa(  1, 3, 5)  -  Qa( 0, 2, 4)  (18) 

fio  =  -  (/o  (1  -  /1)  (1  -  h)  h  (i  -  Si)  (i  -  h))  - 

fo  (1  —  /1)  fl  (1  —  S3)  Si  (1  —  fb)  + 

(i  -  fo)  (i  -  /1)  h  (1  -  S3)  (1  -  Si)  Sb  + 

(1  -  fo)  fl  (1  -  h)  f3  (1  -  fi)  fb 

J  =  circ[(— 1  +  f)2  /2,  —  ((— 1  +  f)2  /3)  ,  (—1  +  f)2  f2  (—1  +  2  /) , 

(1  -  2  /)  (-1  +  f)2  f,  (-1  +  f)2  f,  -  ((-1  +  /)2  /2)  ] 

A0  =  0 

Ai  =  0 

A2  =  2  (1  +  e)  (-1  +  /)3/2 

A3  =  -6  (-1  +  /)2/3 
A4  =  2(l  +  e)2(-l  +  /)3/2 

A  5  =  0 

There  are  only  three  zero  eigenvalues,  so  this  odd  chirality  model  is  sufficient. 

3.3  FHP 

na  =  ^Qa(l, 4)  +  ^Qa( 2, 5)  -  Qo(0, 3)  +  Qa(  1, 3, 5)  -  Qa( 0, 2, 4)  (19) 

fio  =  -  (fo  (1  -  Si)  (1  -  S2)  S3  (1  -  fi)  (1  -  fb))  + 

(1  -  fo)  fl  (1  -  h)  (1  -  h)  U  (1  -  fb) 

2 

fo  (1  -  /1)  f2  (1  -  f3)  fi  (1  -  fb)  + 

(1  -  fo)  (1  -  A)  /2  (1  -  h)  (1  -  /4)  h 
2 

(1  -  f0)  fi  (1  -  f2)  f3  (1  -  fi)  fb 
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2  ;2 


J  =  circ[(— 1  +  /)  / 


((-1  +  /)2  f2  (1  +  /))  (_!  +  ff  p  (_j  +  3  /) 


(1-2/)  (-1  +  /)2/2 


(-1  +  /)2/2  (-1  +  3/)  -((-!  +  /)2/2 


A0  =  0 

Ai  =  0 

A2  =  (l  +  e)2(-l  +  /)3/2 
A3  =  -6(-l  +  /)2/3 

A4  =  (1  +  e)  (1  +  2  e)  (— 1  +  /)3  p 

A5  =  0 

There  are  only  three  zero  eigenvalues,  so  as  is  well  known  the  FHP  model  is 
sufficient. 

3.4  FHP  Without  3-Body  Collisions 

fio  =  ^Qo(l,4)  +  ^Q„(2,5)-g„(0,3)  (20) 


Oo  =  -  (/o  (1  -  h)  (1  -  h)  h  (1  -  h)  (1  -  h))  + 
(1  -  jo)  h  (i  -  h)  (i  -  h)  h  (i  -  h) 

2 

(1  -  /o)  (i  -  h)  h  (1  -  h)  (1  -  h)  h 


J  =  circ[(l  -  f)  f 


.2  (-1  +  /)3/2  (-1  +  fYf 


3  f2 


,(1  -f)  f 


3,2  (-l  +  /f/2  (-1  +  /)3/2 


Ao  =  0 

Ai  =  0 

A2  =  (l  +  e)2(-l  +  /)3/2 

a3  =  0 

A4  =  (1  +  e)  (1  +  2  e)  (— 1  +  /)3/2 

A5  =  0 

There  are  four  zero  eigenvalues,  one  to  many,  so  as  is  well  known  the  FHP  model 
without  3-body  collisions  is  insufficient. 
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4  Conclusion 


A  simpler  model  than  the  original  two-dimensional  FHP  lattice-gas  automaton 
for  simulating  subsonic,  viscous  hydrodynamics  does  exist.  It  is  interesting  that 
although  the  chirality  of  its  collisions  is  definite,  no  spurious  invariants  appear  in 
the  macroscopic  limit  and,  furthermore,  that  this  lattice-gas  automaton  is  time 
reversible  invariant.  The  observation  that  a  fixed  chirality  lattice-gas  model 
can  have  the  correct  macroscopic  limit  is  interesting  for  the  two-dimensional 
hexagonal  lattice.  However,  in  the  case  of  the  three-dimensional  face-centered 
hypercubic  (fchc)  lattice  [3],  this  observation  should  have  more  practical  value 
in  helping  to  reduce  the  number  of  collisions.  The  are  24  nearest  neighbors  in 
fchc,  so  a  full  collision  table  has  224  entries.  For  parallel  computers  such  as 
the  CAM-8  and  CM-5,  implementing  such  a  large  collision  table  is  inefficient. 
Compression  of  the  fchc  collision  table  has  been  explored[4].  Fixing  the  fchc 
collision  table  chirality  should  allow  additional  compression. 

I  would  like  to  thank  Dr  Bruce  Boghosian  of  Thinking  Machines  Corporation 
for  a  related  discussion. 
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